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ABSTRACT 


A nonlinear digital model of a "D" type marine boiler is 
developed with emphasis on water level shrink and swell 
phenomena. Results indicate that further specifications 
of the phenomena are necessary before operational values of 


drum level are attained. 
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IMEUE WBRODUCILION 


Recent developments in control theory and practice 
coupled with digital electronics advances have produced 
increased interest in digital power plant modeling. Although 
conventional marine boiler modeling is not a new concept, it 
is still an area of increased interest due to the relative 
Er analytical data. This deficiency is due primarily 
to the extreme environment that must be endured by sampling 
equipment. Advances in data collecting and processing equip- 
ment have made boiler data collection practically feasible, 
and in fact, nuclear steam generator development depends 
heavily on actual data. However, it is not economically 
feasible for a single manufacturer to retrofit modern data 
collection systems on a standard marine boiler. The majority 
of marine boiler models available are not attractive from 
the control engineer's viewpoint because of one or more of 
the following: 

a. The model is over-Simplified to the point of being a 

"teapot" model. 

b. The model is completely developed and written in 
Laplace transforms and the state space equations are 
too difficult to extract. 

с. The model is over-complicated to the point of being 


computationally inefficient. 





d. The causes and modeling of "Shrink and swell" 

phenomena are ignored. 

This model = to compromise between simplification 
and complication while including a theory and model for shrink 
and swell. The model follows the path specified by Fini [1] 
restated below: 

a. A general D-type marine boiler model is prepared in 

Continuous System Modeling Program (CSMP) language. 

b. А general FORTRAN program is used to prepare the 

initial conditions needed by the CSMP model. 

rhe Initial condition program depends only on data 

E xewupeusned from the manufacturer's technical 


manual and engineering handbooks. 





II. MODEL CONSIDERATIONS 


O BOILER TYPE 

The D-type boiler was chosen for modeling because of 
its comparative standardization among manufacturers coupled 
with its popularity in naval ship power plants. This 
particular boiler is one installed in the LHA-1 class of 
United States naval ships. Since D-type boilers operate 
under a wide range of geometry and output conditions the 
initial condition and model programs are designed to operate 


over similar ranges. 


B. MODEL NOMENCLATURE 
The model nomenclature differs from reference [1] because 
of a need to clarify the model equations and afford readability 
to the actual computer programs. The following objectives 
were considered when developing the nomenclature. 
a. The equations must be easily read, necessitating 
minimum variable-definition referrals by the reader. 
b. The notation used in the model development must be the 
Same as that used in the computer programs. 
The resulting notation combines a. and b. and consists 
of frequently used notations for variables followed by a 
EMbseript to indicate state point location and/or time using 


Figure 1. 


Table 1, а tabulation of general rules and examples, 
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ЕО suffice to allow the reader to follow the model 


development. 


С. EQUATION DEVELOPMENT FORMAT 

The dynamic model and initial condition programs are 
developed simultaneously with corresponding equations 
written in appropriate proximity to each other. A notation 
ЖЕ С! following an equation indicates that it is an initial 
condition equation, while a notation of [ME] indicates a 
model equation. In some instances the equation applies to 
Men programs, an indication of [IC, ME] is used to indicate 
ENS. 

In as many instances as possible equations are developed 
directly from elementary heat transfer and fluid dynamic 
principles. This stems from a desire to encourage more re- 
search in this particular field and to foster a better under- 
Standing of marine boiler dynamics among students of marine 
engineering. The model was developed to be read and used by 


students with a minimum of additional research. 


D. ASSUMPTIONS 

Various assumptions are made to facilitate either model 
EEUuDIÍICation or program tractability. The latter are not 
Ava ys desirable and further research or different programming 
techniques may obviate the need for them. The following 
major assumptions are made; 

l. The energy transferred to the generating tubes is 


distributed uniformly, i.e. "even heating" exists. 


n 





Be Тпеге is a thorough mixing of the fluid in the steam 
drum. 
3. Steam generation occurs only in the boiling section 
of the screen and main bank risers. 
4. All generating tubes and downcomers in a particular 
circuit have the average geometry. 
Uniform energy transfer is assumed to facilitate model 
simplification and is justifiable since the generating 
banks are designed to optimize uniform heat transfer. Mixing 
of the fluid in the drum is a result of the turbulence in the 
w The third assumption results from a lack of reliable 
experimental data and the knowledge that D-type boilers are 
natural circulation boilers relying on a density difference 
in the circulation loop to promote flow. Steam production 
outside of the generating banks would decrease the density 
Emerence, thus inhibiting natural circulation. The geometry 


assumption is necessitated by model simplification. 





III. MODEL DEVELOPMENT 


A. GENERAL 

The boiler is developed from control volumes for which 
the appropriate equations for energy transfer and storage 
rates are written. The initial condition equations are 
generated by setting the rate (d/dt) terms equal to zero and 
solving for the desired variables. 


The model is developed with three different elements: 


l. Furnace-side heat transfer 
2.  Water-side heat transfer 
3. Water-side circulation 


B. FURNACE-SIDE HEAT TRANSFER 

The gas flow path can be followed on Figure 1 where it is 
represented by the dashed line. The fuel and air mixture 
enters the furnace at state point Q where radiation heat 
transfer occurs to the screen riser tube metal. The flue gas 
then leaves the furnace and exchanges heat via convection with 
the superheater (state point R-S), main bank riser tubes 
(state points S-T), and economizer (state points T-U) in that 
order. 

At state point Q the fuel-to-air ratio and fuel and air 
temperatures are relatively constant. This allows the 
sensible heat of the fuel and air to be lumped. The mass flow 


rate into the furnace is the sum of the air flow and fuel flow. 
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MQQ = MFQQ t MAGG [me | 


The initial values of fuel and air flow at a specified 
operating point are known. Therefore, the corresponding 


initial condition equation is: 


MQQg =MFaag +magag = [rc] 


The mass flow rates through the remaining sections of the 


boiler are the same. 


MOQ = MRR = MSS = MTT = muu [mE] 
M = = = 
баб MRRÓ 552 =MTTZ =Muug [Ee] 
The total heat supplied to the boiler can be written as 


the product of the mass flow rate of fuel and the "lumped" 


fuel heating value. 
QQ = MFQR Жж Ену [me] 


At a specific operating point the heat released to the furnace 
is available from the technical manual as is the furnace 
volume. The total heat input at steady state is the product 


Of these two. 


ав = Сто ж Бием. \re | 


This allows the computation of the fuel heating value. 


FRV = деб Р MFQQS \re | 
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d. Furnace to Screen Risers 
The heat transfer to the furnace screen tubes occurs 


Semmerily through radiation, i.e. 


9 7 СА (Tea aS) 


which is the Stefan-Boltzmann lawlwhere oc is the Stefan- 


Boltzmann constant defined as: 


N = 2.9567 хло "ти 
- Pr? RR 


Sec 


Therefore the heat transfer to the furnace screen can be 


written 


Qi = Sie MAA ((TRR чод. 5) ж ж Ч.Я 
a (туу it: Heg. d) KK 4.65) [me | 


where 


SIG6MAA = COXA 


At a given operating point, technical manual values are 
Specified for the furnace heat absorption rate and the area 

of the radiant heat absorbing surface of the furnace screen. 
The product of these two yields the steady state heat transfer 


rate which in turn can be used to compute SIGMAA. 





i > 

This equation assumes that all the radiation from the combus- 
tion flame strikes the screen tubes, and both the flame and 
the tubes behave as black surfaces. 





QiÓ — QFuRAG Ж ARAD [xe] 


Sigman = Q41¢ / ((ткАр +Чор. Ф) ж ж Ч.Ф – 
(тмүф + чор.д) жж Ч.д) 


An energy balance equation for the mass of flue gases in the 


furnace region can be written in the form 


a (m Co E. total energ Supplied E energy 
ат 


transferre to Screen eS 
Energy of gas leaving пасе, 
5 region 


а (massae. x CAR.K TRR)=GQ-QL-MRR 
dt Ж све. ж (тее-тТАмв) 


Note: TRR-TAMB is the absolute temperature. At steady state 

the rate terms are equal to zero. Values of total heat 
supplied to the boiler, heat transfer rate to the screen 

tubes, total flue gas flow rate, and flue gas temperature 

are available from previous equations. This allows the specific 


heat of the flue gas in the furnace region to be computed. 
све. = (вай - Q19) / (MRRd x (тағф-тама)) Tre] 


During a transient, the density of the flue gas in the furnace 
and the Specific heat of the furnace flue gas can be regarded 
as constants near a steady state operating point. The dynamic 


equation can be written as 
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OTRR = (aq -Q1- M&&. X све. ж (кк. - тамв)) / 
(MASSAR X CAR) [me] 


where 


MASSQR = RFuUUE X FurvoL ucl 


Thus, 


L 
2, 
TRE, - | OTRE = TRRA 
t 


or in CSMP-III language 


тее. = IWTERL (TRAP я Отв.) [me] 


The value for TRRO is available from the technical manual. 
2. Flue Gas to Superheater 

For the remaining tube banks heat transfer occurs 

primarily through convection. The energy given up by the 

flue gas in flowing over the tube surfaces is represented 


by the general equation 


energy — mass flow * specific 
given up rate otf heat of 
flue gas flue gas 


* change in temperature of the gas 


я => m Co AT 
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At steady state, the energy transferred to the superheater 
tubes from the gas must equal the energy transferred from 


the tubes to the steam. 
QU - ачФ re] 


The steady state values for superheater inlet and outlet 
steam temperature and pressure are specified by the technical 
manual along with the superheated steam mass flow rate. This 


allows the direct computation of 040 апа 030. 


GM - MMMÓ X(uuud -uMMg) [ic] 
HMMO is the enthalpy of saturated steam corresponding to drum 
pressure, computed with curve fitted equations [2].  HNNO is 
available directly from the technical manual. The specific 
heat of the flue gas in the superheater region can now be 


computed. 


085 - 36 /(мяғф ж(текф -т550У)) Ee 


Thus the dynamic equation for heat transferred from the flue 


gas to the tube metal is written: 


Q3 =MRR XK ORS 7K (TR&.—Tss ) [me] 


The heat transferred to the tubes via convection may also be 
represented by the Grimson correlation [3] which states that 


for cross flow over tubes; 


ЕБ 





Since h can also be written 


в - a /AAT 


then 


a = АС. шәй)” С s ar (4) 
Ve а 


also 
e 
Us = Mm. 
leading to the general equation 


= AC I ^ p. ls Y ae ar (Ke) 
Уре А 
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tna a particular steady state operating point the physical 
properties can be regarded as a constant and the physical 


dimensions are constant, i.e. 


AC Ра E N Р- 3 = constant 
Ур E А 


which implies 


а = (constant) 4 m" AT 


For flow over tube banks with differing rows of tubes and 
differing tube geometries the constant (n) varies little and 
the average is approximately .6 [3]. Therefore, in the case 
of the superheater and main bank riser tubes, the correlation 


is used in the form 


q= 0% AT 
or for the superheater 
Q3 = KRS X MRR X X.G ж (TRS—TWW) 
where TRS is the average flue gas temperature 


TRS = (Ter +Tss) /2.¢ 


Since the steady state values for 930, MRRO, TRRO, TSSO, and 
TWWO are given by the technical manual the constant KRS can 


be calculated. 


KRS И Q 3d / (RR. dx s „O X% (тез - TW WÓ)) се 


” 
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To provide the dynamic solution of superheater flue 
gas outlet temperature, the two transient equations for the 
flue gas to superheater tube metal heat transfer are equated 


resulting in 


IRS-TWW_ = месі ж CRS /KRS 
тее. – т55 


where 


ТЕ — (vef. « vss /2. d 
Solving the above for TSS gives 


т55 = (тее ж (Рн. - уља, э vw) / (енгіз) [me] 


where 


PHIL = 2, Жж мк жож HR ces/KRS me] 


3. Flue Gas to Main Bank Risers 
The equations for the flue gas to tube metal energy 
transfer in the main bank are derived Similarly. For the 


equation involving flue gas specific heat; 
QS = mss x cesT 3x (TSS "rf [me] 
The Grimson correlation equation for the main bank risers is 


Q5- KST X mSS s XK BO ж.(т<т-туч) 





where 
Tst = (Tss ЕР 


By equating these two equations as was done with the super- 
heater the flue gas temperature leaving the main bank can 


be written 


TTT = (тэх ж (рнга,-л) +2, ж туу) /(енгг,+1) [ме] 


where 

Рнг2, = 2. Ж Mss X ox @.Ч ж сат /кзт [mE] 
The steady state heat transfer is computed differently than 
that for the superheater. Since the steady state heat trans- 
fer to the screen bank has been computed, the energy 
transferred to the main bank is the difference between the 
total energy transferred and the desuperheater and screen 
bank energy transfer. At steady state the mass flow rate of 
steam out of the boiler is equal to the mass flow rate of 
feed in the boiler and the total energy transferred is the 
mass flow rate multiplied by the enthalpy change between the 
outgoing steam and incoming feedwater. Thus, the steady 


State energy transfer equation for the main bank is 


сос = MM ж (нммб = H8S6 )—G2$ -6446 үте 


Q990 is computed in the water side heat transfer section. 
Since Q50 = Q60 at steady state the heat transfer 

coefficient KST and the specific heat CST can be computed. 

Кт = 959 / (MSSD KK D.C K (TST -туур)) Ac] 

CST =as¢ / (м55%5 ж (Tss¢ -TrTg)) uc] 
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4. Flue Gas to Economizer 
The economizer tubes differ from the generating and 
Superheater tubes in that they are finned. A reasonable 
finned tube heat transfer correlation was developed by 


Weirmann, et al. [4] 


J = l Pr #З 
Ce & Mar 





G = maximum mass flow rate 
max 
ds correlation constant 
Eo = specific heat of flue gas 
Since h = AY the above 
AT 


correlation can be written 


9 = 3 ААТ lp Gmax 
Pp 2/3 


En x is a function of mass flow rate into the tube bank and 
tube bank geometry. Similar to the case of the Grimson coeffi- 
cient used for the superheater and main bank tubes the physical 


and geometric properties are considered constant leading to 


the use of the correlation in the form 


Q = Cm AT 
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Therefore, for the economizer the equation is 


ОТ = KTU * MTx Gummer — TAX) 


where 
TTU = TTT Тал 
ur 


Paralleling the superheater and main bank, the heat transfer 


may also be written 


QT= MTT X CTH * (TTT- том.) [mE] 
The steady state equation development follows that of the 


superheater and main bank sections and is not repeated here. 


The equations are: 


QD =MAAD ЖҖ (Hed - HAAG ) Nec: | 
0790 = QIQ [re] 


CTU = QTd / (MrT$ «(rrv$ — Tuud)) [те] 
KTru. 7 Q7$/ (wrr$ X(rrud — TAA Ø) lre| 


TTu¢ = (ттт SER /2. d [re] 


Knowing the values for KTU and CTU developed by the initial 
condition program allows the solution of the flue gas outlet 


temperature in the transient analysis. Thus, 


D) Tu = (TTT * (PHES -1.9)+2.B K Txx)/ 
(Ригз 44.0) 
[ME] PHIS -2.d X @TU/KTW 


24 





C. WATER SIDE HEAT TRANSFER 
1.  Risers 

The heat supplied to the riser banks is used both in 
the boiling and nonboiling length of the tube; however the 
ASSumption is made here that the heat transfer is all under 
fully developed flow boiling conditions. Since the nonboiling 
length of the risers is relatively small, little error is 
produced. For fully developed nucleate flow boiling, Levy 


[5] suggests the equation 


4/3 а 
Q= A Psat Alsat Bee 
1.7% х \0° sec- ft 


By lumping the constant area and denominator the riser equations 


can be written in the forms 


Q = С, ж PSAT "a A NTSAT 


Therefore, the two riser equations are 


G2, - ку о<Атжж (4.Ф/з.0)ж(гу/-тзһтуж «.з.@ [mE] 


бо-КУуж “зат жж(4.9/3.0)ж(туу -төхтуж жз.0 ME] 


for the screen and main bank respectively. A simple energy 


balance on the tube metal yields the metal temperatures, i.e., 


mass X Cpe X AT = B 
ат = 4.- 4.4 
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The equations for screen and main bank metal temperatures are: 


orvv = (q1-Q2) / (mAssvV ж сеу) [mel 
Түу = імтвел (туу, бту/) IME | 
отуу = CQs-@c)/ (Massy Ж СРУ) [me] 
ТУУ - іктвві. (туу, 0тУУУ [me] 


The heat transfer coefficient KY can be computed 
directly from tube geometrical data available in the 


technical manual; 


КУ - Аевлм8 /).1%2.Е. 20 tre] 


Since Q10 = Q20 at steady state, the dynamic equation for 


Q20 can be solved in a steady state situation for KV, 


[Ic] KV=Q1¢ f (Pss X (4/3) (rw -TSAT) Xx 3-9) 


Q60 was computed in Section III-B.3. This allows 
the calculation of the initial main bank riser tube metal 


temperature, 


[1c] TW= (%@/(ку ж Р$АТ жж (Ч/з)) ж Ж (1%) +Т5кт 


In the case of the screen risers, the initial tube 
metal temperature is available from the technical manual and 
the initial heat transfer rate has been previously calculated. 
This allows the calculation of the screen riser steamside 


heat transfer coefficient, 


псі kv = Q4D/(Psav xx 4.64.9) x (rw -TSAT)¥x 3.9) 


26 





2. Desuperheater 
пае Dittus=Boelter correlation [6] is used to compute 
the heat transfer rate from the steam to the desuperheater 


Mies. This correlation in its basic format is 


Nu = .Ø23 ЕТ” 


where 
Nu => Nusselt number 
Re => Reynolds number 
Pr => Prandtl number 


Bes en or cooling fluid) 


When the appropriate variables and constants are substituted 
for the dimensionless numbers, the Dittus-Boelter correlation 


eam be written 


А 


where the constant K is defined as;: 


<= 2,3 ж Thermal conduchvity ж heat krarster area 
tube diameter 


С 
( \ M T Hb | 
\ < Те di ameter XK Viscosity x NDO es 


ж bo 


and AT is the log-mean-temperature difference (LMTD). 


Thus, for the desuperheater, the equation is 
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QA = KUP K% MON KK .% ж LWTONP 


KNP can be computed directly. Thus, 
KNP = ((.Ø23 X THCONN) /ONP) X (u. d / (PE. DNP 


X VISCON X NT&805S)) x X d$. 8 X FRAN) 3€K Q3 
Жж AREADS Lre] 


The energy given up by the steam to the tubes may 
also be written in terms of steam specific heat, flow rate, 


and temperature difference. 
QA = CNP K MNN X (TNN -Tee) ime] 


Similarly since the heat transfer to the water in 


ШЕ (Птпп іс totally by convection; 


649 - Kx (Tzz - uH) me | 


и icon darttron calculations THH is considered 


to be equal to TSAT. Therefore 
ка = 8996 / (таб -тЗАТ) [xe] 


The desuperheater outlet temperature is solved by 
equating the two dynamic equations for Q9 and solving for 


the desired outlet temperature. 


TPP = (TNN ~- T£ Zy/ (exe (xe / (CP MNR Ф.2))) 
FTER [me | 
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The instantaneous desuperheater tube metal tempera- 
ture is obtained by an energy balance on the tube metal 


similar to that for the riser tube metal. 


mass «СожаГ-а - 
ағ МЕ dot 


Solving the appropriate equation for dT/dt results in 


OT2% = (Q4-Q99) /(masse K CPZ) [me] 
TEZ = тек. (ТЕР. оТЕЕ іме| 


The specific heat of the steam in the desuperheater 
(CNP) may be calculated using the known steady state values 
for desuperheated steam flow rate and the desuperheater inlet 


and outlet temperatures and pressures. 


04% - мымб ж (HUNG - HPed y Tec] 
CNP = @64Ф/(мык@ ж (тыма -тею@уу [го 1 


With the value of Q90 computed the initial log-mean- 


temperature difference can be calculated. 
LMTONP = 998 /(kuP Xmuud x x g.v) 


The log-mean-temperature difference is a function of 
steam temperature in, tube metal temperature, and steam 
temperature out; therefore the initial tube metal temperature 


may be calculated. 
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Тф = (тм - exeoos x reed) / (4.6 - Execos [m] 
EXPOOS = ЕхР((Тммб -тбей)/ымтомрРУ) rc] 


3. Economizer 


Paralleling the desuperheater development the Dittus- 


Boelter correlation is used to relate the heat transfer from 


the tube metal to the feedwater, 


ШІ (725 NER 


TRUS, 


the only difference being that n is now .4 vice .3 because the 


fluid is now being heated. The appropriate constants are 


again lumped yielding 
QV - KX X MAA X хк 0.9 Ж LIMTOAG 


An energy balance on the water gives the heat absorbed in two 


other forms; 


QS = MAA x% (ugeg - НАА) 
0% = МАА Ж CAR ж (тев — TAA) (ШЕЛ 


1 


At the specified operating point the inlet and outlet condi- 
tions of the economizer as well as the mass flow rate are 


given, allowing the computation of the steady state heat 


transfer. 


Qed = mand x (Head -Haad) = L te 
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The heat transfer coefficient KX is computed directly. 


кх = ((.623 KTHCONA) /OAB) * (4.¢ /(PLX DAB XK 


VISCOA * NTUBEC)) KX Ø. T  PRAA CX di. 
XK AREAEC [re] 


This permits the solution of the steady state log-mean-tem- 


perature difference. 


LMTOAR = ©к@ /(кх ж MAAGD X X "NS [s | 


Following the identical path delineated for the desuperheater, 


the tube metal initial temperature can be computed. 
Txxd - (TAA - Tas x EXFOEL)/(1. 8 -EXPOEC) Erel 
ЕХРОЕс, = ЕхР ((1686 — ТААФ ) /ьмтодв Ice] 


An energy balance on the economizer tube metal yields 


its instantaneous temperature. 
pTxa = (q2-Q7)/(massx ж оғх) [me] 
TAX = INTGRL (TAG, DTXX) [mel 


By equating the formulations for Q8 involving LMTD 


and specific heat, the economizer outlet temperature may be 


calculated. 
TAB = (rAA -Txx)/ (exe (xx | (GA x MAAK 9.2 ))) 
+ TXX [me | 
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The specific heat of the feedwater in the economizer 
is calculated at steady state using the computed value of 


Steady state heat transfer and the given mass flow rates and 


inlet and outlet temperatures. 


САВ = aed /(mand x (тз8@- ТАА@)) ire. 


4. Superheater 
The development of the superheater equations follows 
that of the desuperheater and economizer and for that reason 


the development will not be repeated. The equations are 


listed below. 


быз (мы Жаммм ж (тым-тым) [mel 
ade E MMMO ж (Hung - Нммд) (ro | 


CMN <а80 /(мммф ж (TNNd-Tmmg)) Lee] 


The initial superheater log-mean-temperature difference 
NUxconmputed directly allowing the subsequent calculation of 


the tube-metal-to-steam heat transfer coefficient used in the 


Dittus-Boelter equation. 


Lm Tom = (TNNd-TmmB)/ (Aloe ((Twwg-Tmmß) / 
(Y wwó - TNNÓ ))) [re] 


KW = Q4d /(mmmg x %0.% жамтомы) Wel 
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Ап energy balance on the tube metal yields the instantaneous 
tube metal temperature, the initial tube metal temperature 


being available from the technical manual. 


Drww = (83-04) /(massw xx ору) Сме] 
TWW = INT6RL (Twwg,0TwWw) іме] 


The superheater outlet temperature can now be calculated. 


TNN 7 (rm —T ww) / (Exe (kw /(oma K MN ж 
$.2.))) HF TWW [me] 


D. WATER-SIDE CIRCULATION 
l. General 

The water/steam side circulation equations are by far 
the most difficult to visualize and codify. These equations 
must be accurate in order to predict phenomena such as shrink 
and swell while appropriate assumptions and simplifications 
must be made in order to make the equations tractable. 

As can be seen in Figure 1, the feedwater enters the 
economizer at state point A, passes through the economizer 
to the steam drum where it becomes part of the water volume. 
The liquid leaves the drum via the downcomers, both main bank 
and screen at state points G and C respectively. The main 
bank downcomer delivers its liquid to the main bank risers 
Via the water drum, and the main bank risers deliver the fluid 
to the steam separators. The flow through the screen risers 
is the same with the exception that there is no water drum in 


this circuit. The steam separators separate the majority of 
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the water from the steam leaving the drum at state point M 

and the majority of the steam from the water being discharged 
back into the drum liquid. The steam, with a negligible 
amount of water passes through the superheater via state 

point M-N. At the outlet of the superheater the steam leaving 
at path III is used for "main" steam; that leaving via path 

II travels to the desuperheater where it is cooled, then 
leaves via state point P and is used for auxiliary steam. 

The water leaving the separators with a small amount 
of entrained vapor is mixed with the incoming feed water and 
ШІ а “Тосатіпр"! vapor/liquid mass in the bottom half of 
the drum. The liquid from this mass leaves via the downcomers 
emis Circulated through the loop. 

2. Downcomer Pressure Drop 

A closer look at the downcomer flow is necessary to 
justify assumptions that are made in the development of the 
pressure-drop equations. 

Circulation ratio 1s defined as the ratio of the total 
weight of steam liberated to the drum [7]. For a 600 pound 
poc borler this circulation ratio is on the order of 20:1 
EN Ud 3]. This implies that for every 21 pounds of liquid 
flowing down the downcomers, 20 pounds of it has already 
traveled up the riser and is at Saturation temperature. There- 
fore, assuming the downcomers are perfectly insulated, it is 
reasonable to assume that the downcomer liquid is at or very 


near saturation temperature. 





The momentum equation for the downcomers may be 
written: 
pressure at top - pressure at bottom 
= frictional loss - gravitational head 
+ entrance loss + bend loss 
+ exit loss + inertia force 
The inertia force term may be considered negligible [9]. 
This is the quasistatic approximation which basically states 
that pressure waves move much more rapidiy through the system 
than the important time constants. This is mathematically 
equivalent to the elimination of a large negative eigenvalue. 
This quasistatic approximation is only good for on-line 
transients and does not apply for extremely large discon- 
tinuities which would result in the boiler being taken out of 
ED on (e.g., a ruptured tube). Therefore the downcomer 
momentum equation can be written: 
pressure at top - pressure at bottom 
= (friction factor * downcomer length 
+ downcomer diameter + entrance factor 
репа factor нете factor) 

* (downcomer mass flow ES (2 * (dowvncomer 
cross sectional ned * density of fluid 
in downcomer * gc) - density of fluid in 
downcomer * gravitational acceleration * 


height of downcomer + E 
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In the model notation 


Pepo = (FO X 1L0/0c0 + ENTRCO + BENOCO 
+ EXITCO) ж мее, x 2.0 / (2.6 KADK KAD 
ж &нсер ж GC)— RHoCO 3€ & *K zCD /&C. 


papa = (есєн  169/06H +ENTRSH + GENDGH 
+ EXITGH) Ж м6 XCK 2.5 / (2. X ASHXX ZØ 
Ж Кноєн ж GC) -Еновн ж 6 ж z&n /ac 


for the screen and main bank downcomers respectively. FCD and 


КОШ ате friction factors of the form 


Ба а / (4.24 - 2 106 (R/ks)) 


where R is the pipe radius and KS is the relative sand rough- 


ness [11]. 


FCD 


1.6 /( 4.14 -2.6 kx aos ió (ksco) [reme] 
4.d [ (4.nd -2. 6 Ж ос (ksqH)) Lre,me [ 


FGH 


3. Riser Pressure Drop 
The momentum equation for the riser boiling section 
must take two phase flow effects into account because the 
flow in the boiling section is not homogeneous. Thus, there is 
a relative velocity between the liquid and vapor phases here. 


The steam separators are included in the riser length. However, 
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the effective length of the separators and thus the pressure 
drop are considered negligible because of a general lack of 
information concerning these items. 

The pressure drop due to single phase or homogeneous 


flow in the nonboiling riser section can be written 
= : T Е т | 
AP SDF AP acceleration AP friction AP ravity 


o & 


MPacceleration zum ж (e a а on) 
ас ж А" Sor Kin 


қатан 





J 
E 
ж 
+> 
x 
9” 
5 
+ 
Ж 
Se 
N 





Qe DX AA 


АР асау; = Ss Ж. Һе ж (е. ДЕ Е | 
2, Ж Qc 


Experiments conducted by Babcock and Wilcox [10] 
indicate that these homogeneous flow pressure drops may be 
modified to give the appropriate two phase flow pressure drops 
memo correction factors that are functions of outlet quality 
and operating pressure. The two phase flow pressure drop can 


be written; 


E: * 
AP opp M cereo acceleration 
+ * : Я 
ск "friction 

AP * тг. ) 
gravity gravity 


where the r terms are two phase flow correction factors 


available by using curve fits of the data from reference [10]. 
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Ehe form of these г terms become: 


RGRAVE. = 2:1. 194 xe  жж2.6-0.5000 ж ХЕЕФ +.977% 1] 
Р6РАМК. =2Ч.194 Ж хы. X X 2.6 -05.5 к х2 +.%ТТ© [rc] 


RAcLE - IS.4S CX XFRÍ X x 2.6 dai] x x FF -. 6607 [Exe] 
RACUT = 1974564 Ж Xu x ed 18.494 X Xd —. ddéd7 Tre] 


RFRICE -—-34. 0822, % ХЕРД Ж Ж Р.Ф +23. 7164 * ХЕЕФ 
ee 59734 е. 


RFRICK = -34. 6922, * XUB KK 2.2 +23, 104 K XUD [rel 
RE 


Therefore, in the boiling region of the riser 
= 2 & л a 
A P+ pF = N m. + Ж + (boiling length m Пас 
4. S, А Ac (diam ) £4 А = 
= Che: ht of Den region) Rin p rev 
Ae 


The total pressure drop across the length of the riser is; 





AP = AP oor + АРтрғ 


In model notation the equations are 
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ғорЕ> (моожж2.Ф ж(бноев-енсоо))/ (60, Ж RHCEE. ж 
Кнооо Ж АСЕ. ж ж 2.4) --(Ч.6 Ж FDE. X ООЁ 
Ж моо к ж 2.5) / (&c K DDE x (Rico + 
RHOEE) ж ADE Ж Ж 2. 5) + (< жоє ж (Єнсро 
+RHOEE))/ (SC X2.8) + (moo x X 2.6 X 
RACLE.) / (&C X RHOEE. K AER X K 2.) + 
(2.5 % FEF X LEF моо ж Ж 2.6 Ж RFRICE)/ 
(SC Ж DEF KRHOEE K AEF X X 2.5 + 


(6 Ж жЕР Ж RHOEE ж. KSRAVE) / GC 


PSPL = (мтр y Ж2.0 ж (Енокк.-Еноттуу/ ес. Ж RHOKKOK 
RHOTT X ATK Ж 2.2) (H.P K FIKK MIT HH 
29 $1TK)/ (Ge % OTK K (RHOIT + RHOKK) 
АТК ж Ж 2.4 +9 Ж атк ж (ноті +кнокк)/ 
(66 3 2.5) (mar x 2.@ Ж RACLK)/ 
(GC K RHOKK X AKL X X 2.2) + (2.6 Ж 
Ғы. Жік жмтгжж2.6 Ж кныск)/ 
(Gc. X OKL K RHOKK Ж Ақы ж ж 2.4) + 
(6 Ж ZKL Æ еңокк ж RERAVR) / GC 
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Шр ГЛ оп гастоге аге in the same form as those 


for the downcomers 


Ғов- 1.0/(4.4-2.0 Ж дісві6 (кое) 1 
БЕР = 1. " (1.74 -20 Ж Ао619 (к<еғ)) Lre] 
Ктк- 1.0/(4.74-2.0 Жмевід(кӛтеу) (Тед 


FKL = 1.9/(4.74-2.0 Ж А.ов40 (кӛ к1)) Yuc] 


At a specified steady state operating point the 
pressure drop across the downcomers must equal the pressure 
drop across the risers and the downcomer flow rate must equal 
the riser flow rate, Therefore by equating the appropriate 
pressure drop equations and solving the resultant relation for the 


flow rates, the initial flow rates may be computed as 


меғбф-((кнфеод жө -жғсо -6 Ж агефж((Кнегод 
FRHoEEZß)/2.8) —-G Ж ЕЕЕФ Ж КНОЕЕФ 
RGRAVE) /((FCOK LCO/OCD +ENTRCO + 
BENDCO + EXITCO)/ (2.9 KACOHK 26 
Ж RHOCOP) +((RHoEED -RHODD” ) / 
(КновеФ Ж кнсооб ж Арт ж Ж 2.4)) + 
(4.2 № FOE X LOEØ X О гү x DDF 
X (RhoEE +еноро@) Ж AOF Xx 2.5) + 
RACLE /(RHOEED X ADF X Ж 2. dr (4. Я 
Ж FEF X LEF X RFeICE/ (2.6 :« DOF 


A RHOEED CK. ADE OK 3 2.0 y) xc 2.5 [re] 
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маф = ((RHOCHD KEK ARGH - GK ZTKPSX ((RHOTTU 
*tRHOKKÓ /2.2) - G ж E KL K RHOKKA K 
RGRAVK ) /((F&H x оён /Овн + ENTESH + 
GENO&H 4 Exir&H)/ (2. d X A&H X X 2.0 
Ж. КНН +-(єнокк@ - еногг@ )/ 
(Кнокк@ ж оюнолдё жт. жож 2 @)) + 
(4.9 ж ғтк ж Kg X 2.5)/Q. 6x 0TL 
Ж (енокк@ + кноз:гй) ж Ат. ж ж 2.@) + 
RACLT / (RHokKd KX ATL X X 2.0) -4- (4.d 
& FKL X ж веску / (2.5 жогк 
Ж Кноккф ж Ат. жж 2.6))) X Xd.s rel 


eee Hiser Continuity 
The relationship between the riser inlet and outlet 
mass flow rates is written in terms of the continuity 


оон for one dimensional unsteady flow, 


= oVeA = - а (оа 
ах > At 5 


eV 

In model notation, this becomes 

MLL = MIT- DRHOSLX VOLTL [me | 
MFE = MOD — DRHODF X VOLDE [me] 


The numerical differentiation technique used by CSMP- 
ІІІ іс highly suspect, as are other techniques. Ап "averaging" 
system is used in the actual dynamic model. In addition, the 
flow rate down the downcomer and into the riser is held con- 


stant for the open loop boiler model. 
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o. Riser Quality and Density 
The average density in the risers must be solved for 

КООШО ЛУ Ту. Linearly varying quality along the riser tube 
length follows directly from the assumption of uniform heat 
flux along the riser tube length. The average density along 
the tube length is the sum of the total change in density and 
the density entering divided by the riser tube length [1]. 
Assuming the density varies only in the boiling length of the 


riser, the total change in density can be written; 


.e(4)44 


beilina length 


The average density is 


8 average iD 8 T Fe % поло hing ergh 
В, 


where - density of saturated water 


"f 
L = total tube length 


 шпсе p(l) varies linearly from entering to exiting, p(l) can 


be written in the form 


i 
есе B Ve + Хоз CF UNG) Veg Lt 


Le 


specific volume entering 


where Ve 


Pow 


quality at riser outlet 
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L, - boiling length 
LN, = nonboiling length 


V = change in Specific volume from saturated 
Пао а со) saturated Vapor 


met tegralot the explicit equation for p(2) may 


be solved after rearrangement to 


(V; = S [ме Vig үр, Xo. Vig N 
Us Lg 
ШЕ іс of the form l 


e 


This yields the average density, 


EN Le m e Uta i 
Sav - к Vta | Ve 


ar Ene | 
EN OOCI notation the formula is written 
RHODFO® = (1.8 /Lor )* (verd /((xrFg ) X Vá) 


* Агов (((хғғФу/уе)жуға ж4.0) 
T &uo DOÓ :K LDEÓ) Үге,ме| 


&wen d. - (A. d / cO) (ud / (Occo) лға) 
Alos (( Xu? жуға) /уе +^.6) + RHOJID 
* Uk Ute ‚me | 
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To permit calculation of the quality term the energy 


КООШО С оп the liquid in the riser tubes is evaluated. Thus, 


rate of change = flow energy - flow energy 
of riser energy in out 


+ thermal energy in 
a (ehv) = Min hin 7 Mout Mot TES 


If the enthalpy term is treated as the average enthalpy in 


the riser the equation may be written 
Fo Chg thot /АУУ Арьен 


model notation this is written in the form: 


Онхы. = (мев ж(нтт- не - ха. ж нға /2.4) + Q@ -muL 


Ж Хі Ж HFA /2.ф\ [ (RHOTL X момка \ [me | 
HALL = INTERL (нхо , OXLL) Еме] 
yut = (нх. – пБ) ж 2. ф/ ғ ме 


OuxEE = (мее, Ж (Ноо НЕ -ХЕЕ Ж не /2.6) +92 


-МЕР Ж ХЕБ Ж не6/2.6)/(Енооғ A VOLEF) [mel 
HXFF = INTEEL (HAFFZ,OHXFF) [mel 


ХЕР = (HXFE - НЕ) ЖЗ. 6/n FE [me] 


The initial conditions, HXLLO and HXFFO are calculated in the 


initial condition program. 


The steady state equations for quality are derived 
iterently. The energy transfer rate to the riser liquid 


may be written in the form 


e = mAh 
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where 


Ah = h out =h., 
= he + Хо} hta -- E 


Esdvons for Sr yields 


Xout = С + (he hin )) / (5 hte } 


ХЕЕЯ = (81$ + МЕЕФ № (ВЕ -Н009)) /(меЕФ к нев) [xc] 
KLLD =(QSP + MULE җ (ие - HTT) / CLUB k HFa) [го] 


6. Riser Boiling Boundary Location 
Escunonbori:nge length of riser tube is the product 


of the riser tube length and sensible to total heat ratio [1] 


The sensible to total heat ratio may be expressed as the 
enthalpy change in raising the water to saturated conditions 


ШЕ ес by the total enthalpy change, i.e. 


ЖЕ тан he Б hen 
de hig Xot hga ) — hin 


The equations for nonboiling length are 


LOE -— LOF *K Amax4 ((ur -uoog), 9.8) /(( 4008 


-хЕРб Жиға)-нооФ) tre] 
LIKD = LTL K AmMAXL (CHE - HTIY) B98 )/(( use 
TXLLd % нға)- нта) їгс] 
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The Fortran function AMAX] is used here because the 
above equations are used in an iterative loop in the initial 
condition program and there is a possibility during iterations 
Of reaching a situation where the enthalpy entering is greater 
m UNSabturation enthalpy. The final initial condition solu- 


tion prevents this. The AMAX] function is not used for the 


dynamic model. Thus, 
LOE = LOF X (HF -H00) /((HF +XFF ж HFá)- Hoo) me] 
LIK = LIL X (HF - ux3)/(C HF ха ж не) HIT) [me] 


The boiling volume is derived directly from the solution for 


the nonboiling length. 


boiling length 


boiling volume = total volume * o Б Непсе, 
VOLEF = VoLDF x LEF/LDF [me | 
VOLKL = VOLTL Ж LKL/TL [me] 


7. Steam Drum Liquid Mass and Energy Balance 
The rate of change of liquid mass within the drum is 
equal to the sum of the mass flow rates entering and leaving 
the drum. The liquid mass flowing into the drum is the 
Saturated liquid from the risers, the liquid from steam con- 
densation in the drum, and the incoming enter: The 109019 
leaving is that leaving via the downcomers. The dynamic model 


equation is 
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DOMASL = MLL ж (^.d - xc.) + MEF X (4.d - xrr) 
+MCONO 4- M8S — cc - M66 [me] 


The instantaneous mass of liquid in the drum is 


DMASL = INTERL (oOmasıg, Doması) [mE] 


where 
OmASLD = (уы.оем Ж RobpemL)/2.d re] 


Bir Е ОЕ НЕ balance is derived similarly i.e., 


Ормони = мо ж (^.d - x.) НЕ МЕРЕ Ж (1.6 -xrr) 
* POL — MEG XK HEE Imel 


The instantaneous drum energy is 


DM ORL = INTERL (omong DMOHL) [^E | 


hae initial drum energy is the product of the initial 


imam Liquid mass and initial drum enthalpy. 
Omon Ld = DMASLG ж HDRumg ire] 
Пеи 1 а1 атит Liquid enthalpy is 
рн = Отон. Готаз 


Іп the dynamic model, drum liquid enthalpy is considered the 


enthalpy of saturated liquid. 
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The steam condensation rate equation is based on the 
difference between the pressure and temperature of the steam 


ШІ (Шат оғ the liquid [12] 


DONO = 500.43 E (Emm ¡Gram ANDES x 2.5 


—PSAT / (TSAT - Ho. D )-K X Ф. 5) 02568 nd 


The water level in a marine boiler is generally con- 
sidered in reference to a mid-drum zero, that is, a water 
level of plus one inch implies the water level is one inch 
above the drum centerline. For computational efficiency the 
assumption is made that for small changes in water level 
around the midpoint of the drum the water surface area remains 


constant. This allows a simplified level equation, i.e., 


LEVEL = (Omov - VOLOR M / 2.9 ) И (LSsTMOM, ж ostmom) [me] 


DMOV is the equation for the total volume of "liquid" in the 
drum. Recalling that a small percentage of steam leaving the 
separators is entrained in the liquid, the rate equation for 


"liquid" volume in the drum is 


Dpmov = (( MEF + MLL-MGG ) Ж УЕ +ред Ж (мы. кмге) 


Ж уу 4 "88 x ves —(mec +m64 ) k VA 
Trato x уға) [mel 
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where DMDVO is the volume occupied at time zero which is half 


the drum volume. 


EECUrculation System Initial Condition Iteration 
Procedure 


The initial conditions for the circulation system 
are found by flow rate balancing in the downcomer-riser flow 
loops. An initial guess of flow rate is made using the assumed 
riser exit quality of .05 percent which is reasonable for this 
type of boiler [7]. Coupling this assumption with the assump- 
Anat the percent carryunder is zero, an initial approxi- 
mation of the flow rate can be determined using the figures 
for initial heat transfer rate to the risers calculated in 


menon III-A. Hence, 


heat transfer rate 
assumed quality * latent heat of 
vaporization 


mass flow rate = 


A first approximation of downcomer enthalpy for both 
banks can be calculated. This downcomer enthalpy is assumed 
to be the same for all downcomers. USing an energy balance 


on the liquid in the steam drum 


МЕЕ@ = Q1¢/(XASumE KHFG) [rel 
E = Q5¢/ (XASumE K HFS) Vc] 


HCOD ((MFFd *mu.ó - need x ue + mead ¥ Heed) / 
(MF Fo м LLE) [re] 


ненф = Heo” Lr] 


И 


The screen riser inlet enthalpy is assumed to be equal to 
that of the screen downcomer; however the main bank fluid 
absorbs additional energy from the desuperheater located in 


that circuit, Therefore, the main bank riser inlet enthalpy 
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must be calculated separately. Thus, 
нту@ = НН - Gad /munug үте] 


The downcomer density must be calculated for use in the pres- 


sure drop calculations. Hence, 

Увоф = ((мРЕФ 3-мі.0-мазафуж уе - wea x 6ag)/ 
(mtr mug) Lie] 

RHotog = 1.0 /veod се| 

RHOGHD - KHoCod ro] 


The riser outlet quality is calculated using the 
initial quality formulation previously developed. This com- 
putation is followed by the calculation of the riser non- 
КООШ се апда boiling lengths. 

The average density of the risers is calculated along 
with the two phase flow multiplication factors for use in the 
та 2tLe/pressure drop calculation, 

An updated flow rate is now computed and compared 
With the first approximation. If it is within a specified 
Seger criteria calculation stops. If not the previous 
approximation is updated and the calculations continued with 
the new approximation. 

Upon completion of the flow rate balancing the initial 
flow rates in the downcomer/riser loops are known along with 
the riser outlet quality. 

The steady state drum specific volume and density 


е 


are computed as 
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VORMLP = 0002 + род ж yy 
Коркть = 1.9 /VORMLE 


ihe initial steam mass in the drum is then calculated as 
05тмФ = VOLORM X Сноу /2.6 


9. Superheater Pressure Drop 
All entrance, head, and exit losses are considered 
megligible in the superheater compared to the frictional 


pressure drop. 


AP = fL mf 
A" Qaverage 


By lumping the constants the pressure drop equation is 


PmM=PNN = KNL K mmm ЖЖ 2.7 / Киомьу 


Тһе following definitions apply here: 
Kon 3 = (Pmmø + PNG) /(RHOMME +RHONOA) uel 
КНоммф - (Еномм @ + RHONNG ) /2-Ø 

Kon =(Cemmd – еммф уж Riomug ) /MMMP * K 2 [тс 


10. Steam Valve Equation 
The flow through the steam valve is considered 


EBrectly proportional to the outlet pressure and valve 


Epening, i.e., 
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Ф 
-— 


m= CXPX (valve Opening) 
ммм Е = РМА) Ж коюў Xx vave \me] 


NUEEconstant is computed in the initial condition program. 


KONH = mmm TII / Nh. ved *K PANG) MON] 


E. EQUATIONS OF STATE 


The equations of state listed below are used in both the 
initial condition and dynamic programs. With the exception 
of the subcooled specific volume equation used for the 
feedwater entering the drum, they are reasonably accurate 
in the 300-1500 psi range. The subcooled specific volume 
equation is accurate in the 6ОО-%СО pse range. 

PSAT = ExP ((Alos (Hsar)-4.40788)/.20452 Yel] 
TSAT = ЕхР((.22151 ж ALoa (PSAT) +4.977123)) (21 
HSAT = ExP (P.2G4S2 A ALOE (Psa) + 4.467¢7) 122 
HFG = 922.15 – @.09<510 Ж SAT tA тт -ФЧ 

k PSAT KK 2.6—-4.219E-¢e KK з. ln] 
КЮЕ - (63.6 — d. D1TEA 3& TSAT р 1.132.Е -ØS 


Ж Т5АТ X x 2.5 — G.'TeGE - x K TSAT 
Ж ж з.д \2 | 


von = .ф!еўрч%%# -.@#Дд И Дх@\4ф ж твё 
ғ 3.05ҺЕ-0% Жж твв ж Ж 7.@ 
- $. 1ЧХЕ-П Ж 166 %Ж 3.6 +1-ЧРЯЕ-1З 
X TBB K ЖчЧ.ф - 1. 143Е-\ © K TEE XX SP 
t 9.034€ -29 X T8G CX, р 
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БҮР RESULTS AND CONCLUSIONS 


A. GENERAL 


 іесіпсе of the initial condition program and the dynamic 
boiler model program is given in Appendices A and B respec- 
Bee. The initial condition program output data must be 
properly formatted for input to the CSMP-III dynamic model. 
In addition, because the dynamic model utilizes the liquid in 
the steam drum as a saturation state point from which all 
other state points are derived, two of the initial conditions 
must be modified to eliminate a discontinuity between the 
Steady state program and the dynamic model. The equations 


involved and an explanation are given below. 


DM ASLD = (voLorm ж, &нсоб@ ) /2-@ ucl 
DMOHLH = DMASLÍÉ ж HORUMD (ге. 
DMDHL = INTERL (DMDHLZ , DOMOHL ) [me] 
DMASL = INTERL (DMASLD ,DOMASL) [me] 
DH = омон /OMASL [mel 


As stated previously, the enthalpy of the liquid in the 
drum is considered saturation enthalpy. The dynamic model 
must begin with the rate of condensation in the drum (MCOND ) 
as close to zero as possible. This is a natural steady state 
Pesiation— drum liquid and drum steam both at the same pressure 
and temperature. In order to insure this to be the case 


DMASLO and/or DMDHLO must be modified such that the initial 
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drum enthalpy (DH) very closely approximates the enthalpy of 
ENuUrESted liquid corresponding to the drum Steam temperature 
and pressure (PMM and TMM). This is easily facilitated by 
the use of the CALL DEBUG statement in the dynamic program 
With this procedure: 
l. The initial condition program is executed and the 
initial conditions formatted for CSMP-III use. 
2. the dynamic model is executed for only a short гип 
time, i.e. 5 seconds. 


Observing the DEBUG output from the model, DMASLO 


C2 


and/or DMDHLO are modified such that PSAT and TSAT 
equal PMM and TMM. 

The model is reexecuted for a short run to check, 
MCOND should be very small. 

9. Since DH is implicitly related to PMM and TMM, the 
procedure may have to be repeated a few times. The 
objective is to force as many of the "DERIVATIVE" 
terms in the DEBUG output to comparatively small 
figures as possible. 

The model as used operates satisfactorily with DMASLO - 
geez. 2 and DMDHLO = 3.96163E 06. This forces the initial 
MCOND term to .35949. 

Because of the CSMP function DERIV used in the program, 
the model is extremely sensitive to integration time step. 
Numerical differentiation is not a desirable function to 
perform in a dynamic model; however attempts to explicitly 


differentiate the equations concerned failed. One solution 





Ene problem is to "smooth" the derivative function by 
averaging it over several 22 The model performed 
satisfactorily with a fixed-step integration procedure 
(Runge-Kutta), an integration time step of .04 seconds, 


and averaging the derivative over sixteen time steps. 


B. OPEN LOOP RESPONSE 

The open loop response to a ten percent increase in 
throttle opening is shown in Figures 2-7. AS expected, the 
response of the main bank circulation loop and screen bank 
 СИТагіоп loop is different. This is in keeping with the 
different purposes of those two loops, steam generation and 
mimemeace Screening respectively. The effect of the valve 
Opening increase is barely noticeable in the screen circula- 
pun loop. 

The swell effect is not noticeable. Further conversa- 

tions with Mr. Paul Weitzel at Babcock and Wilcox indicate 
пе percentage of "carry under" is not à constant as It 
was treated in this program. One to two percent is а good 
starting estimate for steady state; however during a transient 
"carry under" mass flow rate is computed by subtracting the amount 
of steam leaving the boiler from the amount of steam produced. 
This should always be a positive number. The vapor that does 
not leave the boiler is "carried under." The following program 


Eus implement this "carry under": 
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Modify 

the section of the dynamic 
model titled "COMPUTE THE 
DRUM SPECIFIC VOLUME" to 
read: 


PROCEDURE CRYUND=FILTR8 (MMM ) 
IE TIME GT.0.0)GO TO 45 
CRYUND=PCU* (MLL+MFF ) 


GO TO 46 
45 CRYUND=(MLL*XLL+MFF*XFF ) -MMM 
46 RISEl-DELAY (250,RISTIM,CRYUND) 


RISE2=CRYUND 
ПЕСМА ШОЕ ОТО 51260 TO 55 


RISE=RISE2 
GO TO 57 
55 RISE-RISEI 
57 CONTINUE 
ENDPROCEDURE 


DDMDV-( (MFF*MLL-MBB)*VF-CRYUND*VV... 
*MBBQ*VBB-(MCC*-MGG)*VF-MCOND... 
*VFG-RISE*VV) 


Follow the above statements with the unmodified 
equations for DMDV. 


The previous procedure was not implemented in the present 
dynamic model. The quality formulation in the present model 
is apparently too simplified and this causes a shortfall in 
quality at the outlet of the screen riser bank. As a result of 
this shortfall, there are instances during transient operation 
when the model is not "producing" as much steam as it is "using". 
Mr. Weitzel suggests using twenty node finite difference 
approximation for quality. This could be done using an equation 
mere quality such as: 

X 2. 9. 2 

A Se is hts 


where 


vá 
І 


quality 


distance up the riser (ft) 


N 
| 


o6 





MESE input го riser (BHUJ/s-ft) 

A = cross sectional area of riser (ES 
оғ = density of saturated liquid ( 1bm/ft?) 
E. ое воет. of liquid entering- (t/s) 
E parent heat of vaporization (BTU/lbm) 
Simultaneously an implicit equation for continuity could be 


applied to produce a balanced mass flow rate without assuming 


that the downcomer mass flow rate remains constant. 


fee ONCLUSIONS 

The model presented is not a finished model. Further 
research is needed to successfully implement the shrink and 
swell theories presented. When the previously mentioned 
сосу with riser outlet quality is solved, the model 
should approximate the dynamics of a wide range of D-type 
marine boilers depending on the initial conditions supplied 
ШО е initial condition program. | 

The initial condition program develops the initial condi- 
tions necessary for the detailed boiler model with the 
relatively scant data available from the boiler technical 
manual shown in Appendix C and a small amount of data from 
common engineering handbooks. Because of the relatively 
small amount of data required „пели та л сопатстопй program 
and/or the model can serve as a basis for a boiler condition 


monitoring system. 
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D. SUGGESTIONS FOR FURTHER RESEARCH 

It is evident that the model suffers by using the CSMP- 
NE tunetion DERIV. In order to facilitate the explicit 
solution of the continuity equation for the riser banks some 
form of differentiation is required, either numerical or 
explicit. An explicit solution is much more desirable from 
a Stability standpoint. А more detailed investigation with 
fewer assumptions might produce the correct form for the 
explicit differentiation equation. It should be noted, 
however, that when using explicit differentiation the model 
must average the derivative over at least two time steps to 
avoid the creation of an algebraic loop. This is easily 
done with the same PROCEDURE format used to average the 
derivative in the current model. 

The use of small perturbation techniques to linearize the 
model and thus allow a state-space representation for multi- 
variable control development and analysis should be undertaken. 
Mere is a computer code available locally to facilitate this 
for CSMP models; however it does not accept the CSMP DERIV 
function and/or the averaging procedure. A means of bypassing 
this problem would result in a general D-type boiler model - 
applicable to a very wide range of boilers currently in use - 
which could be linearized about specific operating points. 
These individual models could then be used for the development 
and testing of multivariable control systems. 

Using locally available optimizing computer codes, an 
optimal D-type boiler design could be attempted with regards 


to boiler geometry. 
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A more detailed analytical and experimental investigation 
of current D-type marine boilers should be undertaken locally. 
Little research has been done in this area and that which has 
been done has often been based on either an incorrect physical 
model of a marine boiler or on a nuclear steam generating 
plant. A starting point could be the data analysis of common 
horizontal and vertical steam separators followed by an opti- 
mum design for these elements. The improper and/or mainten- 
ance of drum internals apparently grossly affects boiler 


operation and water level stability during rapid transients. 
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Figure 4. 
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TvV=INTSRL(TVVC+DTVV 

TYY=INTGEL (TYYU,CTYY) 

т 

ж ENXMPOUTE THE HECAT TRANSFER FOR THE D ESUPERHEATER 

ж ESDUESV T) TUSE FETAL 

* 

09 =С М2 *ЧММ + ( ТЧМ-ТРР ) 

ж . 

t СОУРСОТЕ THE TEMFERATURE OF THE STEAM LEAV ING THE 
2 CESUPEP HEATER 

IPs (TNN-TZZ ) /(E XP (KP /LCN? 98NN50.2)1 e TZ22 

ж SE1 THE 2JOWNCOJNER EN RANCE AND EXIT TEMPERATURES 
* EJUAL TO THc SATUSATI ON TEMPERATURE CORRE SPONDIAG 
> TC CRUN ENTEFALPY 

TGG=TSAT 

THH=TS AT 

TC C=TSAT 

JOC=TS AT 

ж COMPUTE THE ENERGY TRANSFER FOR THE DESUPERHEATER 
f РЕ METAL TC DRUM LIQUID 

BP? KZ*L T22-Tbh) 

* БӘУЕПТЕ ІНЕ RATE EQUATION FOR DESUFERHEATER TUEE METAL 


TEMPERATURE 
TZ22(09-999) /(L 4A SSZ *CPZ) 

COMPUTE THE DESLPERHEATER TUBE METAL TEMPERATURE 
Z221 NTGRL(7 2205 CT 22) 


CO^PUTE THE EZONO'"*IZER ENERGY TRANSFER-TUBE METAL 
NO FEEL WATER 


82MAZASCRES(TEE- TAA) - 


ШЕІРІСГЕ THE FEED TEMPERATURE AT CUTLET 
CF ECONOMIZER 


8Bs(TAA-TXX)/(EXP(KX/ (CAO $448 €x C,2) )) e TXX 


Mele ТӨС RATE cOUATI CN FOR THE ECONOMIZER TUEE 
METAL TEMPERATURE 


TXX2(08-071/( 4ASSX*CP x) 
COMPUTE THE ECONCMIZER TUEE “ETAL TEMPERATURE 
X3sINTS RL (TXXC,CTXX) 
CCMFJT7TE THE SUPERHEATEC STEA^ OUTLET TEMPERATLRE 
NNs (T4M -T4W )J/( EXP (K&W/ (CAM 34NN 50.2) 2) + TWh 
ССАРЦТЕ THE SUPERHEATER ENERGY TRANSFER-TUBRE METAL 
4 = Шы. 


COMPUTE THE RATE EQUATICN FIR SUPERHEATER TUBE METAL 
TEMPERATURE 


ThW= (03-04 )/(*ASSNW2C FW ) 
COMPUTE THE SUPERHEATER TUGE METAL TEMPERATURE 


$ REPO RARA ss eee OH HE ta a X n9 ЧЕЗ HHO He 
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ТМЧчЧ-ЕХР( ,22151%31206(94М)»%8.77123) 
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* ЕКОЕ ХЕ FOR THE ECCNOMI ZER FEED SUTLET EITH4SLPY 
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HB E= (MEB*HALD+C8)/MRB 

% = 

+ FEECRATE EQUATION 

ж 

“BB=MB80 

ж 

ж TWO PHASE FLCwW YULTIOLIER EQUAT IONS 
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RGRAVESZ^.T9^ *XFF**2.)-6.50CÓOÓ*XFF*.91To 

EESANNE24. 719*^*XLL**2.0—€.50€6* XLL*. 9776 

САС| Ес15.456%%ХҒҒ>т2,Әж1о. 4544%ХЕҒ-., С0007 

RACL Jz15.^56-^ *XLL **2,.) tl3 ,*94«* *XLUL—-.00007 

"FRICES-34.Ct22* XFF**52.,C* 23. 7164* XFF «4,8734 

RFRA(K2-73«4.086zcz*X.LL *2Z .CvZz2.Tl164*XLL *.8734 
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ж STATE POINT EQUATIONS 
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nm SAT=DH 

Fr=HSAT E 

ЕКЕс-922.15-0.%2516%Р54741.717Е-0%4%95АТ%%2.С-4.219Е-06... 
*PSAT**3.0 

FVzHS^T *FFC j _ 

kr CF =63 .-8-0.21781%75 АТ +1 2132 E-05 *TSAT #32 -0-6.786Е-02.. o 
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G. Anticipated Performance 


Ra! of Operation - Per Cent 

Toil Steam Generated lb/hr 

Suprheated Steam Ibs/hr. 

De: perheated Steam Ibs/hr. 

Borr Drum Pressure psig 

Surheater Outlet Pressure psig 
Suprheater Outlet Temperature ^F 
Derperheated Steam Outlet Pressure psig 
De:perheated Steam Outlet Temperature °F 
Ect omizer Inlet Pressure psig 

Ect omizer Inlet Temperature °F 

Ect omizer Outlet Pressure psig 

Ec omizer Outlet Temperature °F 

Cang Air Inlet Temperature ^F 

To Air Flow lb/hr 

rot! Oil Flow lb/hr 

Ай ред Efficiency 5» 

ouanteed Efficiency % 

Rantion and Unsccounted for Losses 45 
Exess Air % 

caron Dioxide с 

Naiber of Burners in Operation 

[hittle or Non-Throttle of Air Doors 


Ога Loss - Total Inches Water 
Ihragh Double Casing 

Mhtagh Burner Register 
hragh Boiler & Superheater 
hragh Economizer 


34$ 'emperature Leaving Superheater Screen ^F 

Jas emperature Leaving Superheater °F 

34$ 'emperature Leaviug Main Bank °F 

348 emperature Leaving Economizer °F 

leg Release KB/Hr/Sq. Ft. Radiant Heat Absorbing 

ur ce | 

lea Release KB/Hr/Sq. Ft. Total Heating Surface 

leasRelease K B/Hr/Cu. Ft. Furnace Volume 

"urice 1eat Absorption KB/lIr./Sa. Ft 

lea Absorption First Water Screen Row KB/Hr/Sq. 
Et. (Max.) 

lead Absorption Maximum KB/llr./Sq. Ft 
(Furnace Screen) 


Endurance 


2 Boilers 


56 
195,000 
155,000 
40,000 
690 

665 
912 

651 

635 

715 
250 
703 

355 

100 
242,229 
14,333 
85.7 
85.7 

99 

15.0 
13.0 

6 

N.T. 


13.74 
1:33 
2.60 
4.83 
4.98 


2471 
1826 
694 
208 


447 
15.5 
183.9 
128.4 
210 


210 


APPENDIX C 


Desigii 
Rated 
Full 
Power 


83 
290,000 
235.000 
55,000 
690 

630 

914 
606 

Gog 

726 
250 

705 

B7 

160 
568,092 
21,781 
84.4 
84.4 

99 

150 
13.0 

6 

N 


33.955 
3.09 
6.00 
12722 
12.60 


2594 
1961 
774 
427 


680 
23.6 
2792 
ЕЛІ) 
231 


231] 


Dr Losses (Full Power) Based On 241 Cu. Ft. of 100%F Air/Lb. of Fuel Oil 


*D ft Losses (Overload) Based On 2-41 Cu. Ft. of 100%F Air/Lb. of Fuel Oil 


**laft Losses (Overload) Based on 260 Cu. Ft. of 68°F Air/Lb. of Fuel Oil 


&7 


Maximuin 
Intermittent 


Power 


D 
340,000 
265.000 
75,000 
690 

610 
980 

558 

680 

733 

250 

703 

372 

100 
433,559 
25,656 
Go 
83.8 

9225 

|520) 
15.0 

6 

М.Т. 


47.96 
4.28 
8,50 
18.14 
17.04 


2642 
2021 
812 
453 


800.6 
28.0 

202 
ІРІСІ 
255.3 


2553 


Endurance Boiler 
] Boiler Overload 
114 ІРО 
400,000 420.000 
300,000 320.000 
100,000 100.000 
690 690 
580 570 
898 895 
480 470 
690 690 
749 754 
250 250 
708 710 
3797 380 
100 100 
515.154 540,952 
30,364 32,009 
83.0 82.7 
83.0 82.7 
1.01 98 
15.0 15.0 
13.0 13.0 

6 6 

N.T. NT. 
69.51 76.22** 93 487** 
6.00 10:22 
Е 1.52 
225072515 35.52 
26.40 28.40 29.06 
2686 2699 
2085 2103 
856 872 
454 498 
948 999 
32.9 34.7 
389.6 410.8 
213.19 2215 
288 256 
288 256 





ll. Anticipated Metal Temperature Degree Е. 


Water Screen Tubes Outside 

Water Screen Tubes Inside 

Superheater Tubes Outside - Maximum 

Superheater Tubes Inside - Maximum 

Maximum Inner Casing Temperature 

Maximum Outer Casing Temperature 
at way of Structural ties 


Endurance 
2 Boilers 


Maximum Outer Casing Temperature 


Tube Data 


Side Wall and Roof 
Rear Wall 

Front Wall 

Screen Bank 


Superheater 


Main Bank 
Main Bank 
Economizer 


Desuperheater (in Water Drum) 
Downcomers . 
Downcomers 


Downcomers 


Risers 


Design Maximum 
Rated Intermittent 
Full Power Power 
656 661 
554 554 
1046 1049 
1020 1019 
825 835 
350 350 
145 
O.D. M.W.T. Хо. 
25 134” 71 
2 134” 54 
2 Tds 22 
2 21347 102 
1852 .120 72 
2 .134" 34 
1” 085" 2808 
2 1807 182 
D 220” б 
85/8” 483” 1 
10 3/4” 5197 5 
12 3/4" .601 2 
6” 500” 7 


Endurance Boile: 
] Boiler Overload 
667 680 
554 555 
1048 1045 
1014 1012 
850 870 
350 350 
145 145 
Material 


MIL-T-16286 CL. 
MIL-T-16286 CL. 
MIL-T-16286 CL. 
MIL-T-16286 CL. 


MIL-T-16286 CL. 


MIL-T-16286 CL. 
MIL-T-16286 CL. 


MIL-T-16286 CL. 


5А-268 ТР-430 
Schedule 80 Pipe 
ASME SA-106-D 


Schedule 80 Pipe 
ASME SA-106-B 


Schedule 80 Pipe 
ASME SA-106-B 


ASME SA-106-B 


*Tliese tubes may be MIL-T-16286, Class A (Seamless) or MIL-T-17188 (Sesmed) electric resistauce welded. 
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